Coulomb drag between quantum wires with different electron densities 



Thomas Fuchs and Rochus Klesse 
Universitdt zu Koln, Institut fiir Theoretische Physik, Ziilpicher Str. 77, D-50937 Koln, Germany 



Ady Stem 

Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

(August 5, 2004) 

We study the way back-scattering electron-electron interaction generates Coulomb drag be- 
tween quantum wires with different densities. At low temperature T the system can undergo a 
commensurate-incommensurate transition as the potential difference \W\ between the two wires 
passes a critical value A, and this transition is reffected in a marked change in the dependence of 
drag resistivity on W and T. At high temperature a density difference between the wires suppresses 
Coulomb drag induced by back scattering, and we use the Tomonaga-Luttinger model to study this 
suppression in detail. 



PACS 

I. INTRODUCTION 

A constant electric current in a conductor establishes 
a stationary non-equilibrium distribution of electron mo- 
menta. Relaxation processes with electrons in a second 
("passive") conductor nearby can result in a significant 
drag force, which leads to a measurable drag current 
or drag voltage in the passive conductor. This effect is 
known as Coulomb drag if the momentum transfer is me- 
diated by the Coulomb interaction [1]. The drag effect is 
very sensitive to electronic correlations and can therefore 
be used to probe them. 

While the Coulomb drag between two-dimensional 
electron systems has been extensively studied experimen- 
tally [2], the observation of this effect in strictly one- 
dimensional (ID) systems still presents a challenge and 
experimental work on this subject is quite limited [3-5]. 
The ID case is especially interesting since here electronic 
correlations are much stronger than in two- or three- 
dimensional systems. This is reflected by the fact that 
the Fermi-liquid theory, which applies in three dimen- 
sions and marginally holds for two dimensional systems, 
generally fails in one dimension. Instead, essential fea- 
tures of a correlated ID electronic system are captured 
by the Tomonaga-Luttinger (TL) model [6]. 

A theory for the ID Coulomb drag that is based on 
the TL model [7,8] predicts a behaviour that qualita- 
tively deviates from the one in higher dimensions: At 
temperatures T below a certain energy gap A the drag 
resistivity po between long wires of equal electron den- 
sity increases exponentially with decreasing temperature, 
pu cx expA/T [9]. At temperatures well above the 
gap the drag resistivity shows a distinct power law be- 
haviour, pd on , where the exponent x is determined 
by the forward scattering part of intra- and inter-wire 
interaction. For vanishing forward scattering interac- 
tion X equals unity (and the linear T-dependence of a 
Fermi-liquid approach [10] is reproduced), while it de- 



creases with increasing repulsive electron-electron inter- 
action. The gap energy A corresponds to a correlated 
ground state, where the electrons in the two wires or- 
der in a zigzag formation with period equal to the Fermi 
wavelength Xp. The gap energy A therefore strongly de- 
pends on the backscattering (2kp) component of the ee- 
interaction. It is exponentially suppressed if the distance 
d between the wires is large compared to As- 
under the strictly linear spectrum assumed by the TL 
model a drag response is generated only by inter-wire 
back-scattering. A recent analysis [11] by Pustilnik and 
collaborators addresses drag induced by forward scatter- 
ing in a wire with non-linear dispersion. In particular, it 
is shown in that work that at low temperatures forward 
scattering gives a contribution p^^^ oc to the drag re- 
sistivity if the difference in the Fermi velocities 5vp is 
small compared to T/kp. For 5vp larger than T/kp, 
the contribution due to forward scattering decays more 
rapidly, oc T^. 

This communication focuses on the backscattering 
drag between wires of different electron densities. In prin- 
ciple, the difference in densities is expected to suppress 
inter-wire 2kp scattering, since the momentum 2kp is dif- 
ferent in the two layers, and to affect drag also through 
the difference in Fermi velocities in the two wires. Here, 
we neglect the difference in Fermi velocities between the 
wires (a difference which induces corrections of the or- 
der of ^2YT/i2' '^^^^ l^i electro-chemical potential of 
the i'th wire), and focus on the effect of the difference 
between the relevant 2kp's (a difference which induces a 
much larger correction of the order of ). We envis- 

age the experimental situation of two identical wires on 
different electro-static potentials [3,5]. 

We start by analyzing Coulomb drag in the high tem- 
perature regime, where inter- wire back-scattering may 
be treated perturbatively. In Sec. Ill we reinvestigate the 
drag in the incommensurate phase with different electron 
densities and present detailed results for the non-linear 
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drag at finite temperatures. Parts of these results are 
already published in the review [4]. 

Following that analysis, we address the low temper- 
ature regime. For equal densities, the low temperature 
regime is that of an interlocked crystal, characterized by 
a drag resistance that diverges at zero temperature. As 
an electro-chemical potential difference W is turned on 
between the wires, the relative density between the wires 
is initially incompressible, until an incommensurate- 
commensurate (IC) transition [12,13] takes place [14]. 
This transition markedly influences the drag. The tran- 
sition from the incommensurate to the commensurate 
phase is, roughly speaking, a readjustment of charges into 
an energetically more favorable zigzag configuration. 

We employ the theory of the IC transition [12,13] and 
use a simple Drude model for solitonic charge transport 
to investigate the drag at low temperatures T <C A. We 
conclude that for \W\ < A the drag resistivity is 
exponentially large, pn cx cxp(A — |VK|)/T. Above the 
transition point the drag shows an inverse square root 
behaviour, po oc [W"^ — A^)~^/^, which changes to an 
1/|H^| dependence for large \W\ > A (c.f. Fig. 4). This 
is studied in Sec. IV. 

We hope that the present work may help to better 
understand present and near-future drag experiments. In 
particular, together with the recent results of Pustilnik 
et al. [11] it may help to clarify the role of forward and 
backward scattering in the drag effect. 

II. COULOMB COUPLED DOUBLE WIRE 

We consider electrons in two parallel, strictly ID wires 
of length L separated by a distance d. We assume the 
wires to be identical, and assume that a voltage W/e is 
imposed between them, leading to a difference between 
their clcctro-chcmical potentials = /i2 — pi- It is 
assumed that L is much larger than all other relevant 
lengthscales of the problem, i.e. the wires are practically 
of infinite length. For the beginning we neglect the elec- 
tron spin. 

Using the same notations as in [8] and assuming that 
the electro-static potential is constant along each wire, 
kinetic and potential energy of the electrons are given by 

rwk 

The index r corresponds to left (r = — ) and right (r = -|-) 
moving electrons, w refers to active {w = 1) and passive 
{w = 2) wire. The electron number in wire w is = 
-h A_^, where N^w = J2k ^lw{f'^)0'rw{k). The Fermi 
wave- number kp is defined via the identical equilibrium 
density in both wires for vanishing external potentials. 

The expressions for forward scattering and backward 
scattering electron-electron interactions Hf and Hh, 



resp., are precisely as in Ref. [8] Sec. IIIB, and are there- 
fore not repeated here. 

We proceed as in [8] with a standard transformation 
[6] to bosonic fields 

<l>w{x) = — - V [Q+w{q) + e-w{q)] 

Li Q 

g 

with their conjugated 
q 

where Jw = N^yj — N-y,, and Qw is the electron density in 
wire w. Note that in the present case the particle number 
may change with variation of the external potentials 
and therefore the zero- mode terms must be included here. 
In the following, the inverse momentum cut-off a is set 
to the Fermi wavelength 27: /kp- Then, transformation 
to symmetric and antisymmetric (— ) charge modes 
0c± = 2~-'^/^((^i ± (/)2) decouples the Hamiltonian into 
independent parts, H = Hc+ + He- , with 

^AEo /■ (,Qg _ f dxV^4>c-dx, (1) 

-Ka J 2tt J 

and a similar expression without cos-term for Hc+- The 
parameters u = Uc- and K = Kc- are determined by the 
interwire and intrawire couplings gi and gi as explained 
in [8]. A = gi/2'Ku denotes the dimensionlcss interwire 
backscattering coupling, Eq = u/a \s oi order the Fermi 
energy. 

III. HIGH TEMPERATURE DRAG 

At high temperatures the two wires are only weakly 
correlated, and drag may be calculated perturbatively, 
in a method similar to that employed in [7,8]. Here we 
will focus on the case of wires with different densities. 

In the weakly coupled regime it is appropriate to switch 
to a statistical ensemble with fixed electron numbers 
= Lriy, (w = 1,2). This enables us to define peri- 
odic bosonic fields by 

(t>c- {x) = (j)c- {x) + qx , 

where q = 2-^^^TT{ni - na) = 2-1/2(^^1 - fc^a)- The 
dynamics of ^ follows from Eq. (1) to be given by the 
Hamiltonian 
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H / COS •\/8(0c- — qx)dx + const. (2) 

7ra J 

We calculate the drag using the formalism of Refs. [7] 
and [8], whereby the drag voltage is 
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2kL 



{dxn-)i 



■{dlcl>c-)i. (3) 



Here, n_ = (ni — n2)/2 and n = K/iiru arc density 
and compressibility of the anti-symmetric channel, and 
the averaging is taken with respect to a state carrying a 
current / in the active wire (1). 

Using the condition of stationarity and the equation of 
motion for (j) we find that to second order in the backward 
scattering the drag ££> = eVo/L is 
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(4) 
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Im J dtj dx{cos{\/84>x,t — Qx + Lot) sin(\/8^o,o 



where Q = u5kp, and the term cut = — t results from 
a Galilei transformation of the first wire which generates 
the current / (cf. Refs. [7,8]). 

Calculating the correlator in the usual way, we obtain 
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where the function Nr,^ is given by 



^r,K = ^lirn^ y e"'*''/'' + i) sinhs^ 
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(6) 
(7) 



The result (5) is limited to the weakly coupled regime, 
which requires A ^ T. Further, we have assumed that 
all three relevant energies T, and \Q\ are small com- 
pared to Eq. Boundary effects are also neglected, which 
is allowed as long as the wire length L exceeds the ther- 
mal wavelength ujT . 

As a forward to a discussion of the drag (5) as a fmic- 
tion of the parameters T, Q and cj, wc first analyze the 
function ^r,K- An asymptotic expansion of Eq. (7) 
shows that "^r^K decays exponentially ~ e"'' for large 



positive argument r, and that ^r,K 
negative r. More precisely, we find 



Nr,K w 2 sin(27rii')r(l - 2K) 
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FIG. 1. Nf-./f for interaction parameter K = 1, 0.75, 0.5 
and 0.25. 

For the particular values K = 1, corresponding to non- 
interacting electrons, and K = 1/2 the expression (7) 
simplifies to 



Nr,l/2 



2r 
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While N-r^i bears some resemblance to the Bose distribu- 
tion, Nr,i/2/27r is the Fermi function (cf. Fig. 1). 



A. Non-linecLT drag 

For T < \uj\ and T < \Q\ (but still T > A and 
T ^ u/L ) we can approximate N^-^k in Eq. (5) by the 
asymptotic expression for negative r, and N^^x = for 
positive argument r. This shows that for |a;| — |(5| ::g> T 



sj, = sin2(27ri^)r2(l - 2K) 

neu 



2K-1 



whereas the drag becomes exponentially suppressed for 
\Q\ — ^ T. This result has been already derived 
by Nazarov and Averin [7]. The seeming |a;^ — 
singularity it exhibits near the threshold |cj| = |Q| is 
smeared over a regime where ||w| — |Q|| T, which can 
be quantitatively described by Eq. (5). For large cur- 
rents, \uj\ >• IQI, the drag goes as Ed oc oj^^~^. Fig. 2 
illustrates this behaviour. 



B. Linear drag resistance 



From Eq. (5) one easily finds the linear drag resistivity 

PD = ^\i=o as 
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FIG. 2. Drag in arbitrary units as a function of tj oc / in 
units of IQI = u\SkF \ for interaction parameters K — 0.75 and 
0.25 at temperatures T = 0.25Q (smooth curves) and T <^ Q. 
For large currents (|a;| ^ |<3|) the drag goes as oj''^^^. At 
low temperatures the drag shows an \uj — Q\^^~^Q{oj^ ~ Q^) 
singularity near the threshold ^ \Q\, which is smeared 
over a regime where — \Q\\ ~ T. 




-7.5 -5 -2.5 2.5 5 7.5 

O/T 

FIG. 3. Normalized drag resistivity po (cf. Eq. (9)) as func- 
tion oiQ/T = uSkp/T for interaction parameter K = 1, 0.75, 
0.5 and 0.25. respectively. 



expressions for drag ed and drag resistance pD are ba- 
sically identical to (5) and (8), except of one difference, 
which is essential: The interaction-determined constant 
K is replaced by an effective K = ^{Kg + Kc-), where 
Ks is the interaction parameter of the spin modes (cf. 
[8]). For weak backscattering coupling !v 1, such that 
K = ^{1+K). I.e. the effective parameter K of the spin- 
full system is closer to the non-interacting (Fermi-liquid) 
value Kpi = 1, which reflects the moderating effect of 
the spin-modes. 



(8) 

with po = VSiT^^-^EoX'^/e'^u and N^^^ = dNr^^/dr. 
The drag resistivity is a symmetric function of Q with its 
maximum at Q = 0, where po oc T^^~^. For \Q\ larger 
than temperature the linear drag resistance is exponen- 
tially suppressed, due to the exponential suppression of 
lS5r,K for large and positive r. Thus, the peak of pd{Q) 
has height oc T^^if-s ^j^j jg directly proportional 

to temperature. More details can be extracted from Fig. 
3, where the normalized drag resistance 



~Pd{Q/T) 



Pd{Q,T) 



(9) 



Pd{Q = 0,T) 

is plotted against Q/T iov K = 1, 0.75, 0.5 and 0.25. 

1. Influence of electron spin 

A spin- full system can perturbatively be treated in the 
same way as above the spin-less system. The resulting 



IV. INCOMMENSURATE-COMMENSURATE 
TRANSITION 

Eq. (1) is the Hamiltonian of the sine-Gordon model 

with an additional coupling of the soliton density ng = 
J dx 3x-\/8(/>c-/27rL to the potential difference W. It is 
this coupling which gives rise to the IC transition [12,13]: 
At large AE'o and small \W\ the cos-term dominates 
and suppresses the soliton density to zero, which is 
the commensurate phase. In the opposite situation, the 
cos-term is negligible and the ground state density of 
the solitons is adjusted to rig = ttKW/u. This value of 
rig corresponds to electron densities = (qw) in the 
two wires that are independently tuned by the respective 
electro-statical potentials Vw This is the incommensu- 
rate phase. 

In the classical model, for W less than the energy 
Eg of a single soliton the density ris vanishes strictly 
at zero temperature. The transition to the incommen- 
surate phase sets in at \W\ = Eg, where for small 
\W\ — Eg > the soliton density increases logarithmi- 
cally with rig ~ I ln(|M^| — i^s)|~^, corresponding to a 
critical exponent (3 = [12,13]. At finite temperatures 
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the transition is smeared out, as ususal for ID statistical 
models. 

Quantizing the model modifies the transition in two 

aspects: the transition energy increases to a renormalized 
soliton energy Eg = A and the critical exponent changes 
to ^ = 1/2 [12,13]. 

Obviously, the IC transition is accompanied by a dras- 
tic change in the inter-wire electronic correlations. It 
can be therefore expected that the Coulomb drag varies 
strongly near the transition, and we now address this 
variation. 

In the commensurate phase the electron densities in 

the two wires are inter-locked to one another. The com- 
mensurate phase as a gapped phase has no low-energetic 
excitations and therefore shows no response to a small 
anti-symmetric field Ec- at strictly vanishing tempera- 
ture. Hence, at linear response the resistivity to the flow 
of non-identical currents in the two wires is infinite, and 
so is the drag resistivity. 

At finite but small temperatures T <C A — \W\ soli- 
tons and antisolitons are thermally activated with a den- 
sity n^/a cx exp{-Eg/a/T), where E^/a = A T W' de- 
notes the energy of soliton and antisoliton in the exter- 
nal potential. In the presence of an anti-symmetric field 
Sc- = Si = —£2 both solitons and antisolitons contribiite 
to an anti-symmetric current jc- ■ A detailed calculation 
of this current response to the field Sc-, taking into ac- 
count the damping of the soliton motion is beyond our 
present scope. However, we are able to conjecture that 
the temperature dependence of the anti-symmetric con- 
ductance is exponentially activated, following the tem- 
perature dependence of the soliton density. The conduc- 
tance is the product of the soliton density by their mobil- 
ity. As long as the soliton density is small, we expect the 
mobility to be limited by the interaction of each soliton 
with its environment, rather than by the interaction with 
other solitons. In that limit the mobility will not depend 
on the soliton density, and thus will not be thermally ac- 
tivated. In that limit, then, we expect an exponentially 
large drag resistivity 

PD ~ ae(^-l^l)/^ (10) 

at low temperatures T <C A — |T4^|, where the double wire 
system is in the commensurate phase. The coefficient a 
may still exhibit a significant temperature dependence. 

When the external potential \W\ exceeds the critical 
value A for the incommensurabiltity transition, prolifer- 
ation of non-thermal solitons {W > A) or antisolitons 
(W < — A) sets in with a density increasing as Hg/a oc 
(W^ - A^)^/^ [12,13]. Accordingly, above the transition 
point the conductivity CTc- increases a 
This means that the exponential behaviour of the drag 
resistivity (10) changes to an inverse square root depen- 
dency 




A IWI 
FIG. 4. The inverse drag resistivity of a pair of wires at 

low temperature, T <C A, as a function of the potential dif- 
ference W (schematic). For \W\ less than the energy gap A 
the charges order in a zigzag formation (commensurate phase) 
and the drag is exponentially large, po ~ oexp(A — |M^|)/T 
(I). For W above the gap, \W\ — A ^ T, the drag resistance 
decreases as {\W\ — A)"^''^ with increasing W (II). At large 
W ^ A the drag resistance decays linearly with W (III). 

valid at temperatures T <^ — A. Also the coefficient 
b may exhibit a significant temperature dependence. The 
transition from the exponential behaviour (10) to the in- 
verse square root dependence takes place as W passes 
through an energy window of width of the order of T at 
±A. 

At large external potentials, \W\ A, the 

(anti)soltion density is nearly linear in \W\, such that, 
assuming that the soliton mobility is still independent 
of its density, the drag resistance goes as pd oc 1/|W^|. 
The low-temperature behaviour of the drag resistivity is 
summarized in Fig. 4. 

2. Influence of electron spin 

So far the discussion ignores electron spin and is there- 
fore limited to a spin-polarized system. What will hap- 
pen for unpolarized electrons? Referring to previous 
work [8] the following can be said: The spin fluctuations 
moderate correlations of the charge modes, such that A 
assumes a much lower value than in a spin-polarized sys- 
tem. If r ^ A can nevertheless be realized, we expect 
qualitatively the same characteristics of the drag as in 
the spin-less case. This can be seen by the renormaliza- 
tion group equations of a spin-full system, which suggest 
that both spin and relative charge modes become gapped 
at low temperatures (cf. [8]). Inspection of the backscat- 
tering hamiltonian (see Eqs. (41-43) in [8]) reveals that 
also in this case solitonic excitations in the c— mode are 
possible. This leads to the same phenomenology of the 
low temperature drag as just discussed for the spin-less 
case. 



5 



V. DISCUSSION 

We have seen that the incommensurate-commensurate 
transition of a strongly coupled double wire is ac- 
companied by a marked change of the drag resistance 
(cf. Eq.s (10) and (11)). Observation of this effect re- 
quires that the temperature is of order of or below the 
energy gap A. Taking the estimate of A of order mK 
for a spin-unpolarized double wire with Eq w 10 K [8] we 
conclude that the strongly coupled regime might be elu- 
sive for these systems. However, in case of spin-polarized 
electrons the situation looks much better, where for the 
same parameters one obtains A of order lOOmK . 

The perturbative results derived in Sec. Ill should be 
applicable to existing experimental data [3-5]. In fact, 
both groups observe a peak of the drag resistivity as the 
gate voltages are varied. While the qualitative agreement 
is obvious, a quantitative comparison is difficult because 
the relation between electron densities in the wires and 
gate voltages is not trivial. Eq. (5) suggests that mea- 
surement of the non-linear drag might be an insightful 
probe in future experiments. 
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